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We construct four-dimensional covariant non-linear theories of massive gravity which are ghost- 
free in the decoupling limit to all orders. These theories resum explicitly all the nonlinear terms 
of an effective field theory of massive gravity. We show that away from the decoupling limit the 
Hamiltonian constraint is maintained at least up to and including quartic order in non-linearities, 
hence, excluding the possibility of the Boulware-Deser ghost up to this order. We also show that 
the same remains true to all orders in a similar toy-model. 



Introduction: Whether there exist a consistent ex- 
tension of General Relativity by a mass term is a basic 
question of a classical field theory. A small graviton mass 
could also be of a significant physical interest, notably for 
the cosmological constant problem. 

A ghost-free linear theory of massive spin-2 - the Fierz- 
Pauli (FP) model [l[ - had been notoriously hard to gen- 
eralize to the nonlinear level Q: the Hamiltonian con- 
straint gets lost in general and, as a result, the sixth 
degree of freedom - the Boulware-Deser (BD) ghost - 
emerges as a mode propagating on otherwise physically 
meaningful local backgrounds {e.g., on a background of a 
lump of matter) . Part of this problem can be seen in the 
effective field theory (EFT) approach to massive grav- 
ity [^| in the decoupling limit [3, H| ■ There, the problem 
manifests itself in the Lagrangian for the helicity-0 com- 
ponent of the massive graviton. This Lagrangian generi- 
cally contains nonlinear terms with more than two time 
derivatives. The latter give rise to the sixth degree of 
freedom on local backgrounds, while in general, these 
terms lead to the loss of well-posedness of the Cauchy 
problem for the helicity-0 field theory [J 0] ■ 

A step forward has been made recently in Q where it 
was shown that: (a) the coefficients of the EFT can be 
chosen so that the decoupling limit Lagrangian is ghost- 
free; this involves choosing the "appropriate coefficients" 
order-by-order, and an algorithm was set for this pro- 
cedure to an arbitrary order; (b) once the "appropriate 
coefficients" are chosen in the effective Lagrangian, in 
the decoupling limit only a few terms up to the quartic 
order survive, all the higher order terms vanish identi- 
cally. Moreover, the surviving terms are unique as their 
structure is fixed by symmetries [H, @] . 

In the present work we build on the above two points, 
and go far beyond them. In particular: (1) We construct 
Lagrangians that automatically produce the "appropriate 
coefficients" once expanded in powers of the fields; these 
give rise to theories that are ghost-free automatically to 
all orders in the decoupling limit. (2) Using the obtained 
Lagrangians we study the issue of the BD ghost away 
from the decoupling limit; we show that the Hamilto- 
nian constraint is maintained at least up to and including 



quartic order, hence excluding the possibility of the BD 
ghost up to this order. We also express the exact poten- 
tial for gravity in a simplified (l+l)-dimensional model 
and show explicitly how the constraint is preserved to all 
orders. 

The present framework provides explicit resummation 
of the nonlinear terms in the EFT Lagrangian of massive 
spin-2. Another way to resum these terms is to use an 
auxiliary extra dimension 0, H[. The latter has so far 
been shown to give the ghost-free decoupling limit only 
up to the cubic order In @, H| the resummation is 
obtained via the second order partial non-linear differen- 
tial equation. The present approach achieves this via an 
algebraic non-linear equation. 

Formalism: Define the tensor as the covarianti- 
zation of the metric perturbation, g^ v = rj^ + h^ v = 
+ r] a i ) d l _ l (j) a d u (p b , where the four Stiickelberg fields 4> a 
transform as scalars, and r) a b = (—1,1,1,1), [3]. The 
helicity-0 mode tt of the graviton can be extracted by 
expressing q 



(x a — rf^d^-K), such that 
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d^dyir. (1) 



We may therefore define the following quantity 



K$(g,H) = SS-y/SS-M =-Y d dn( Hn Yu, (2) 

n=l 

(2n)\ 



with 



(1 - 2n)(n!) 2 4 r ' 



(3) 



Here Hj} = g» a H av , and (#")£ = H^H^ ■ ■ ■ H^ 1 de- 
notes the product of n tensors H^. Below, unless stated 
otherwise, all the contractions are made using the metric 
g^ v . The tensor /C M „ = g^aKZ is defined in such a way 
that 



n 



[lis ■ 



(4) 



We use the same notation as in [4{ where square brack- 
ets [. . .] represent the trace of a tensor contracted us- 
ing the Minkowski metric, e.g. [H] = r/^Il^ and 
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[II 2 ] = 77 a/3 77 /il/ n QM II / 3 ; ,, while angle brackets (. . .) rep- 
resent the trace with respect to the physical metric g^ v , 
so that (H) = g^H^ and (H 2 ) = g^g^H^H^. 

We are first interested in the decoupling limit. For 
that, let us define the canonically normalized variables, 
7r = A^7r with A^ = m 2 Mp\ and h^ v = M P ih^. The 
limit is then obtained by taking Mp\ — > oo and m — > 
while keeping 7r, /i M „, and the scale A 3 fixed. First, we 
construct an explicit example of a non-linear theory that 
bears no ghosts in the decoupling limit, and then give 
a general formulation and show the absence of the BD 
ghost beyond the decoupling limit in quartic order. 

Massive Gravity: The consistency of the Fierz-Pauli 
combination relies on the fact that the Lagrangian 



C 



(2) 
dor 



[n] 2 - [n 2 ] . 



(5) 



is a total derivative. To ensure that no ghost appears 

(2) 

in the decoupling limit, it is sufficient to extend C d J T co- 



variantly away from hn 



0. 



replace [II] and [II 2 ] by 



(JC) and (JC 2 ) respectively, so that the total Lagrangian 
reads as 



£ = 



M 2 



r=-g[R-— U {g,H) 



(6) 



with the potential U expressed as an expansion in H as 
U(g,H) = -4((/C) 2 -(/C 2 )) (7) 
= -4(£d n (tf")) 2 -8£d n <#™>. 



(H 3 )) (8) 



n>l n>2 

Expanding this expression to quintic order, 

U(g,H) = {(H 2 ) (H) 2 ) - l -((H)(H 2 ) 

-^((^ 2 > 2 +4( J ff)( J ff 3 )-5(^ 4 )) 

-1 (2(H 2 )(H 3 ) + 5(H) (H*) 7(F 5 » + • ■ • , 

we recover the decoupling limit presented in |5| with the 
special indices c 3 = d 5 = f 7 = 0. 

Note that the Lagrangian © with ([7]) can be obtained 
from the Lagrangian 



M 2 
i-A = — — * 



^(i?-m 2 (/C 2 „-/C 2 )) 



(9) 



where K,^ is an independent tensor field that gets related 
to as in @ due to the constraint enforced by the 
Lagrange multiplier A£. Note, the expression ([2]) can 
be rewritten as = 5% — y/ d fi 4> a d v <j) b r] a i 1 , that gives a 
square root structure in the full Lagrangian. 

Decoupling limit: It is straightforward to notice that 
the leading contribution to the decoupling limit 



-4((n7r) 2 -(d Q cV) 2 ),(10) 



is a total derivative. The resulting interaction La- 
grangian in the decoupling limit is then given by [5J 



with 



M 2 x m 2 5 
8 ShZ 



-gU{g,H)) 



h uv =0 



Using the relations 

SJC(g,H) 

S(IC(g,H) 2 ) 

the expression for X simplifies to 
1 



- (g^ - K,n . 



X 
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(12) 

(13) 
(14) 

(15) 
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The tensor X^ is conserved and gives rise to at most 
second order derivative terms in the equations of motion. 
This tensor can be expressed as the product of two epsilon 
tensors appropriately contracted with powers of 14^,, [6J . 
For the potential (O, the Lagrangian in the decoupling 
limit is then given by, see Ref. [f| 



h^(£h)^ + h^X^ . 



(16) 



and this result is exact (i.e. no higher order corrections). 
Notice that this is also in agreement with the results of [5j 
up to quintic order, for the special case C3 = d§ = f 7 = 0, 
but we explicitly demonstrate here that this result re- 
mains valid to all orders. 

General formulation: As mentioned in at each 
order in the expansion there exists a total derivative con- 
tribution 



4"r(n) = - J2 (-i) m i n [n m ] c (n ~ m) 



(n — m)! 



"'dcr 



(n),(i7) 



with £^.(H) = 1 and £^(II) = [14]. These total deriva- 
tives generalize the "Fierz-Pauli" structure used previ- 
ously to all orders. More generally, the potential of any 
theory of massive gravity with no ghosts in the decou- 
pling limit can be expressed non-linearly as 



U(g,H) 



-4$>„4:)(/C) 

n>2 



(18) 



where [I4 m ] in (|T7| should be replaced by (JC m ) and ex- 
pressed in terms of g and H using Q . 
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Here again this specific structure ensures that the lead- 
ing contribution to the decoupling limit is manifestly a 
total derivative by construction, 



-9U{g,H) 



= total derivative , 



(19) 



and the resulting interaction Lagrangian can be derived 
by noticing the general relation 



8h» 
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(20) 



(21) 



(-l) m n! 



^ 2(r 



ni: 
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c 



(n—m) 



(n). 



using the notation n"„ = r\^ w and = 0. The de- 
coupling limit Lagrangian is then given by (|16j) with the 
same definition (fT2"|) for the tensor X^, giving here 



(22) 



n>2 



with 



where the are the Bernoulli coefficients. We also 
focus on the case where ct2 = 1 and a n = for n > 3. 
In what follows, we work in terms of the ADM variables 

m 



9 



oo 



-TV" 



9o i 



Ni, and = 7y 



(26) 



with the lapse A = 1 
metric 7^ = 5ij 
potential is then of the form 



5N, and the three-dimensional 



hij. In terms of these variables, the 



^gU = A + 5NB + N i N j [-26 ij +C l 



(27) 



-SN(S ij +V ij ) ~ -SN 2 5 l] 



where A,B,C^ and T> % i are functions of hij, at least first 
order in perturbations, and C ij • 2'P" = -\h^ ' +0(h 2 j), 
and in this section we raise and lower the space-like in- 
dices using 5ij. Notice that this is completely consistent 
with the analysis performed in [4, and corresponds to 
setting the coefficients in (43) of[4jto J 4 = £>=Z) = 
E = 0, while C = —1/2. We emphasize here that the 
presence of a term of the form CN 2 N 2 does not signal 
the presence of a ghost, since any quadratic terms in 
the lapse disappear after integration over the shift as we 
prove in what follows. Indeed, in terms of redefined shift 



4" ) = E(-!) 



m=0 



2(n — m) 



-rrm An-m) 
j 11 /^ dcr 



(II). (23) 



This is in complete agreement with the results obtained 
up to quintic order for a 2 = 1, "3 = — 2c 3 , a 4 = —2 2 d 5 
and 0:5 = —2 3 ff. However we emphasize that the results 
in this paper are now valid to all orders. The special 
theory found in 0, S] corresponds to the specific choices 
of coefficients = 1 and a 3 = —1/2, see Ref. [Toj | . 
Furthermore, at each order the tensors Xj^) are given 



by the recursive relation 



x {n) 



(24) 



with xfiS = 1/2^ 



So since Xut = all these tensors 



vanish beyond the quartic one, Xjfiy = for any n > 4, 
and the decoupling limit therefore stops at that order, as 
previously implied in Q. 

Boulware-Deser ghost: The previous argument en- 
sures the absence of ghost in the decoupling limit, but it is 
feasible that the ghost reappears beyond the decoupling 
limit, and is simply suppressed by a mass scale larger 
than A3. Certain arguments have hinted towards the ex- 
istence of a BD ghost, We reanalyze the arguments 
here and show the absence of ghosts within the regime 
studied. To compute the Hamiltonian, we fix unitary 
gauge for which 7r = 0, such that 



(H n ) = Y^{~^YC l i +n ~ 1 \h l+T ' 



(25) 



-5NS& - -5Nti 
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L\m , (28) 



the Hamiltonian is of the form 



K = ^Vl (NR° + NjB?) + (A + B5N) (29) 



m 2 Ml 



Pi T ij 



1 



1 



nuij - -CiUjUk + —n k nirij ) , 



up to quartic order in the metric perturbations. Then, 
it is straightforward to check that the variation of the 
Hamiltonian (|29j) w.r.t. the shift rii gives an equation 
which is independent of N, and serves to determine 
rij. Moreover, the lapse remains a Lagrange multiplier 
even after integration over the shift, hence giving rise 
to a Hamiltonian constraint on the physical variables. 
Whether this constraint gives rise to a secondary con- 
straint, and whether the system should be quantized as 
a first- or second class system, is a separate interesting 
question. The mere existence of the Hamiltonian con- 
straint is sufficient to claim the absence of the BD ghost 
to that order [la ], yet without breaking Lorentz invari- 
ance, fli| . 

The Hamiltonian evaluated on the constraint surface 
is proportional to m 2 and whether or not it is positive 
semi-definite is determined by the explicit expressions for 
A, B,C y and V % K Thus, in general certain backgrounds 
could have slow tachyon-like instabilities, however, this 
is a separate issue from that of the BD ghost that we 
clarified above. 
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(1 + l)-d massive gravity: Proving the absence of 
the BD ghost in complete generality beyond the quartic 
order is a grand task, which we save for a separate study. 
However, we can analyze here a similar issue in a (l+l)-d 
toy-model, where we consider the Hamiltonian 



NR° + 7 n 7Vii?i 



—NU{g,H) 



(30) 



with R and R\ arbitrary functions of the space-like met- 
ric 7xi arl d its conjugate momentum, and the potential 
U is given in (J7J. In 1 + 1 dimensions, it is relatively easy 
to check that the Hamiltonian then takes the exact form 

H = M^^\nR° +~f 11 N 1 R 1 - 2m 2 N\ (31) 



2mf I 1 - •\/(v / 7 + N) 2 - 7 u iVf 



and seemingly includes terms quadratic in the lapse when 
working at quartic order and beyond, 



H~H + HiN + m 2 N 2 N 2 + 



(32) 



By stopping the analysis at this point one would infer 
that the lapse no longer enforces a constraint. However, 
this should be determined after integrating the shift. In 
other words, in terms of the redefined shift ri\ 

jVx = nx (711 + N^) , (33) 

the Hamiltonian takes the much more pleasant form 

H = yfjNR - 2m 2 (1 + ^N) (34) 



(Vl + N)[ n 1 R 1 + 2m^Jl-n( 



which remains linear in the lapse, even after integration 
over the shift. It is again straightforward to see that 
the lapse does enforce a constraint, and does so for an 
"arbitrary background" . 

Outlook: We have given a covariant non-linear real- 
ization of massive gravity in 4D which: (1) is automati- 
cally free of ghosts in the decoupling limit, to all orders 
in non-linearities; (2) keeps the lapse as a Lagrange mul- 
tiplier away from the decoupling limit, at least up to 
quartic order in non-linearities. These findings consti- 
tute what we believe is a very significant step forward, 
and strongly suggests the existence of an entirely ghost- 
free classical theory of massive gravity. However, to prove 
this statement in complete generality, two important in- 
gredients are yet missing: (a) proving that the lapse re- 
mains a Lagrange multiplier to all orders; (b) checking 
whether the secondary constraint is generated or not, and 



whether the theory could be canonically quantized as a 
first or second class system. For the consistency of the 
theory at the quantum loop level one would have to es- 
tablish the existence of a symmetry which protects this 
theory against quantum corrections that could revive the 
ghost. These points will be explored in a further study. 
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The approach of [la ] is equivalent to the EFT approach 
of 0], as was shown in [TJ]- Hence, the claim of [lg] on the 
presence of the BD ghost in the quartic order, if correct, 
would contradict our results. However, what has really 
been diagnosed in [l5[ is the issue already raised in pj, 
which we have just addressed. In particular, the apparent 
ghost-like nonlinear terms identified in [l5| ] , to the extent 
they were presented in [l5| ] , are in fact removable at that 
order by a nonlinear field redefinition, in complete con- 
sistency with our results above. This will be discussed in 
more detail elsewhere. 



